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When a self-gravitating body (e.g., a neutron star or black hole) interacts with an external tidal 
field (e.g., that of a binary companion), the interaction can do work on the body, changing its mass- 
energy. The details of this "tidal heating" are analyzed using the Landau-Lifshitz pseudotensor and 
the local asymptotic rest frame of the body. It is shown that the work done on the body is gauge- 
invariant, while the body-tidal-field interaction energy contained within the body's local asymptotic 
rest frame is gauge dependent. This is analogous to Newtonian theory, where the interaction energy 
is shown to depend on how one localizes gravitational energy, but the work done on the body is 
independent of that localization. 
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I. INTRODUCTION AND SUMMARY 

This is one in a series of papers that develops pertur- 
bative mathematical and physical tools for studying the 
interaction of an isolated gravitating body with a com- 
plicated "external universe" in the slow- motion limit. By 
"slow-motion limit" we mean that the shortest timescale 
r for changes of the body's multipole moments and/or 
changes of the universe's tidal gravitational field is long 
compared to the body's size R: R/t <C 1, where we have 
set the speed of light equal to unity. 

In addition to this slow-motion requirement, we also 
require that the body be isolated from other objects in 
the external universe, in the sense that both the radius 
of curvature 1Z of the external universe in the body's 
vicinity, and the lengthscale C on which the universe's 
curvature changes there, are long compared to the body's 
size: R/K < 1 and R/C < 1. 

The slow-motion, isolated-body formalism, to which 
this paper is a technical addendum, is, in essence, a per- 
turbative expansion in powers of the small parameters 
R/t, R/1Z and R/C. For a detailed discussion of the 
structure of this expansion and its realm of validity, see 
Thorne and Hartle Q. As they discuss at length (their 
Sec. LB), the slow-motion and isolated-body assumptions 
make no reference to the internal gravity of the object un- 
der study. Consequently, the Thorne-Hartle formalism in 
general, and the results of this paper in particular, can 
be applied even to strongly-gravitating bodies, as long as 
the source of the external tidal field is far enough away 
to allow a "buffer zone" where gravity is weak. This 
buffer zone, called the local asymptotic rest frame, will 
be described more fully at the beginning of Sec. III. 

Two examples of isolated, slow-motion bodies are: (i) 
a neutron star or black hole in a compact binary sys- 
tem that spirals inward due to emission of gravitational 
waves; and (ii) the satellite Io, which travels around 
Jupiter in an elliptic orbit and gets heated by Jupiter's 
tidal gravitational field . 

The series of papers that has been developing the per- 
turbative formalism for studying tidal effects in such 
slow-motion, isolated bodies is: 

1. The book Gravitation j|, Section 20.6 (written by 
John Wheeler): laid the conceptual foundations for ana- 
lyzing the motion of such an isolated body through the 
external universe. 

2. Thorne and Hartle Q: formulated the problem of 
analyzing the effects of the external universe's tidal fields 
on such an isolated body, and conceived and initiated 
the development of the perturbative formalism for study- 
ing the influence of the tidal fields on the body's motion 
through the external universe, the precession of its spin 
axis, and its changes of mass-energy. 

3. Thorne m: developed the theory of multipole mo- 
ments of the isolated body in the form used by Thorne 
and Hartle. 

4. Zhang H: developed the theory of multipole mo- 



ments for the external universe's tidal gravitational fields, 
which underlies the work of Thorne and Hartle. 

5. Zhang Q: extended the Thorne-Hartle analysis of 
motion and precession to include higher order moments 
than they considered. 

6. Thorne 0] and Flanagan ||: initiated the study 
of tidally-induced volume changes in the isolated body, 
using the above formalism. Their studies were motivated 
by numerical solutions of Einstein's equations by Wilson, 
Mathews, and Maronetti [|| which seemed to show each 
neutron star in a binary being compressed to the point of 
collapse by gravitational interaction with its companion. 
Thorne and Flanagan found no such effect of the large 
magnitude seen in the numerical solutions. An important 
piece of Thome's analysis came from examining the work 
done on each star by its companion's tidal field — i.e., an 
analysis of "tidal heating." 

Thome's analysis of tidal heating required dealing with 
an issue that Thorne and Hartle had discussed, but 
avoided confronting: For an isolated body with mass 
quadrupole moment Xjk, being squeezed by an exter- 
nal tidal gravitational field Sjk = Rjom (with Rap-yg 
the external Riemann tensor), there appears to be an 
ambiguity in the body's total mass-energy M of order 
5 M ~ TjkSjk- (Here and throughout we use locally 
Cartesian coordinates in the body's local asymptotic rest 
frame; cf. Thorne and Hartle JTJ . Because the coordinates 
are Cartesian, it makes no difference whether tensor in- 
dices are placed up or down.) 

One can understand this apparent ambiguity by exam- 
ining the time-time part of the spacetime metric in the 
body's local asymptotic rest frame: 



2M 2,l ik n j n k 
9oo = -H 1 —s 
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Here n 3 — x 3 /r is the unit radial vector and r is distance 
from the body in its local asymptotic rest frame. Among 
the terms omitted here are those of quadratic and higher 
order in the body's mass M and quadrupole moment Xjk 
and the external tidal field £jk — terms induced by non- 
linearities of the Einstein field equations. Among these 
nonlinear terms is 



^00 



^-jk^jk 



(2) 



whose 1 jr behavior can be deduced by dimensional con- 
siderations. This term has the multipolar structure 
(monopole) /r identical to that of the 2M/r term from 
which one normally reads off the body's total mass, and 
its numerical coefficient is ambiguous corresponding to 
the possibility to move some arbitrary portion of it into 
or out of the 2M/r term. Correspondingly, the body's 
mass is ambiguous by an amount of order 



SM ~ X ik £ 



(3) 
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In Sec. II of this paper we shall see that this ambiguity 
is not a purely relativistic phenomenon. In the Newto- 
nian theory of gravity, there is also an ambiguity of mag- 
nitude (|^) in that portion of the gravitational interaction 
energy of the body and external field which is contained 
within the body's local asymptotic rest frame. In Sec. 
Ill we shall explore the relativistic version of this "tidal- 
quadrupolar" interaction ambiguity by computing how 
the interaction energy changes under a relativistic change 
of gauge (infinitesimal coordinate transformation). 

"Tidal heating" of the isolated body involves inject- 
ing into it an amount of energy that is of just the same 
magnitude as this ambiguity, AM ~ XjkSjk- Does this 
mean that tidal heating is an ill-defined, unphysical con- 
cept? Certainly not, as is attested by photographs of 
volcanic plumes ejected from Jupiter's satellite Io (see, 
e.g., Ref. (!(]]). That volcanism was predicted by Peale, 
Cassen, and Reynolds Q before the Voyager spacecraft 
discovered it; and their explanation — that Io gets tidally 
heated to high internal temperature by the coupling of 
its quadrupole moment to Jupiter's tidal field — is firm 
(see, e.g., Refs. 0JT|). 

In this paper, we use the phrase "tidal heating" to 
mean the net work done by an external tidal field on an 
isolated body. This phrase is slightly misleading, as the 
work done on the body need not necessarily go only into 
heat. The additional energy might be used to deform 
the body (i.e., raise a tide on it) or it might go into 
vibrational energy. If, however, timescales for changes 
of £jk and Xjj, are not close to the body's vibrational 
periods, then, when averaged over many cycles of change 
of £jk and/or Tjk, the work will contribute primarily to 
heat, as in the case of Io. 

In his analysis of binary neutron star systems, Thorne 
argued, but did not prove in detail, that although 
the tidal-quadrupole interaction energy is ambiguous, the 
amount of work done on an isolated body by an external 
tidal field (i.e., the amount of tidal heating) is unambigu- 
ously given by the formula^ 



dW = _1 cEy 
dt 2 ij dt 



(4) 



and he argued that this is true in general relativity the- 
ory as well as in Newtonian theory. In Sec. II we shall 
give a Newtonian demonstration of this claim. In Sec. 
Ill we shall give a relativistic demonstration — showing, 
more specifically, that although the quadrupole/tidal in- 
teraction energy is gauge dependent, the work done is 
gauge invariant and has the value (0). 



II. NEWTONIAN ANALYSIS 

In this section, we consider a Newtonian body, with 
weak internal gravity |$ | <C 1 (where $ is the body's 
gravitational potential), subjected to an external New- 
tonian tidal field. We assume that the external field is 
nearly homogeneous in the vicinity of the body, £ R 
(cf. Fig. 1; in the Newtonian case, the inner boundary 
of the vacuum "local asymptotic rest frame"|] — indicated 
by the inner dashed circle — would be at the surface of the 
body). 

In our analysis, we will consider a variety of contri- 
butions to the total energy inside a sphere which encom- 
passes the body and whose boundary lies within the local 
asymptotic rest frame — i.e., the region where the exter- 
nal field is nearly homogeneous (again, cf. Fig. 1). We 
denote by V the interior of this sphere and by dV its 
boundary. Of greatest interest will be the interaction 
energy (between the body and the external tidal field) 
and the work done by the tidal field on the body. Both 
of these quantities are the result of slow changes of the 
tidal field £jk and the body's quadrupole moment Ijk- 

As a foundation for our analysis, consider a fully iso- 
lated system that includes the body of interest and other 
"companion" bodies, which produce the tidal field £jk 
that the body experiences. For simplicity, assume that 
all the bodies are made of perfect fluid (a restriction that 
can easily be abandoned). Then, for the full system, the 
Newtonian gravitational energy density and energy flux 
can be written as 



H' l '" = ,,(i« a 



n + <i> | + i-*^ 



pv j [ -v 2 + n + - + $ ) - — $ t $ 



(5) 
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where <fr,p,p,v, and II are the Newtonian gravitational 
potential, mass density, pressure, velocity, and specific 
internal energy j?]]. 

Using conservation of rest mass p t + (pv' J )_j = 0, the 
first law of thermodynamics pdH/dt+pv J • = 0, the Euler 
equation for fluids pdv-* /dt+pQj+pj = 0, Newton's field 
equation <&jj — 47rp, and the definition of the comoving 
time derivative d/dt — d/dt + v^d/dx^ , it can be shown 
that Eqs. (||) and (Q) satisfy conservation of energy 



e 



on 
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0. 



(7) 



Equations (||) and (Q) for the energy density and flux, 
however, are not unique. Equally valid are the following 
expressions: 



1 Actually, expression (Q) is only the leading order term 
in the perturbative expansion of dW/dt. The next term is 
— ^BijdSij / dt, where Bij is the "magnetic type" tidal field of 
the external universe and Sij is the body's current quadrupole 
moment || . In this paper we confine attention to the leading- 
order term. 



2 We shall discuss the concept of "local asymptotic rest 
frame" near the beginning of Sec. III. 
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which also satisfy energy conservation (Q) but localize the 
gravitational energy in a different manner from 9^. En- 
ergy conservation (^|) will also be satisfied by any linear 
combination of and O^- Imposing the additional 
condition that, for any acceptable 8, the system's total 
energy 



e = / e 0U rf 3 



(10) 



must be independent of the choice of 00 forces the co- 
efficients to sum to 1. Hence, a perfectly valid form for 
qo m ig 



e ^ 1 = oe?" + (i - a)& 



0/i 
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(11) 



where a is an arbitrary constant. 

Notice that the choice of a gives a specific energy lo- 
calization. For example, a = puts the gravitational 
energy entirely in the field [— (V<i>) 2 /(87r)], so it is non- 
zero outside the matter. This is analogous to the local- 
ization used in electrostatics (1/8-7T times the square of 
the electric field). Choosing a = \ 1 by contrast, puts the 
gravitational energy entirely in the matter (ip$), so it 
vanishes outside the body. When \p^ is integrated over 
the entire system (body plus its companions), the result 
is a widely used way of computing gravitational energy 
(e.g., Sec. 17.1 of Ref. ||). 

The energy in the region V that contains and surrounds 
the body but excludes the companion, 



e 00 d 3 a 



(12) 



will depend on a; i.e., it will depend on where the energy 
is localized. By contrast, the total energy (jnj for the 
fully isolated system (body plus its companions) will be 
independent of a. 

Another way to express this ambiguity of the local- 
ization of the gravitational energy is given by Thornc 
(Appendix of Ref. uy): it is possible to add the diver- 
gence of rjj — f3Q$> } j (where (3 is an arbitrary constant) 
to @1° and the time derivative of —rjj to without 
affecting energy conservation (Q) or the physics of the 
system. Indeed, this method is completely equivalent to 
the one presented above. The constants are related by 
13= (a-l)/47r. 

Throughout the region V, the Newtonian gravitational 
potential can be broken into two parts: the body's self 
field <& and the tidal field <& e produced by the external 
universe (the companion bodies), so that 



$ = $ + $ c 



(13) 



The external field is quadrupolar and source-free in the 
region V so that 



= is 



(14) 



and, furthermore, the tidal field £y evolves slowly with 
time. The body's (external) self field is monopolar and 
quadrupolar and has the body's mass distribution as a 
source: 



M 3Xy? 



$o = 

r 2 



®o,jj = ^P 



(15) 



The quadrupole moment Zjj, like that of the external 
field, evolves slowly with time, but the body's mass M 
is constant. Recall that r = ^/SijX l x : ' is radial distance 
from the body's center of mass and n? = x J /r is the 
radial vector. 

A useful expression for the total energy E\> in the 
spherical region V can be derived by inserting Eqs. (|Tl|), 
(|), (§), and (|13|) into Eq. |L2|). The resulting expression 
can be broken into a sum of three parts — the body's self 
energy E a (which depends only on <I> and p), the ex- 
ternal field energy E c (which depends only on <E> e ), and 
the interaction energy E- m t (which involves products of 
<&o with p or $ ): 



Ey 



E n + E 



E h 



(16) 



where 



E a = 
E c = 
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(17) 
(18) 
(19) 



Inserting Eqs. (|lj) and (|l_5J) into Eq. (|19|) and integrating 
gives the interaction energy inside V as 



_ (2 + a) 



(20) 



which depends on a. In other words, it depends on our 
arbitrary choice of how to localize gravitational energy. 
This is the ambiguity of the interaction energy discussed 
in Sec. I. 

The rate of change of the total energy inside V can be 
expressed in the form 



dEy 

dt 



Q 0j n ] r 2 dn 



(21) 



<)V 



by taking the time derivative of Eq. (|l2|), inserting Eq. 
(p]), and applying the divergence theorem. This expres- 
sion, like that for the energy, can be broken into a sum by 
combining Eqs. (0), (§, (§), (||), ©, ©, and ©: 
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dE v 
dt 



dE c 
dt 



dE u 



dt 



' dt 



+ [ -^[a$ ,t$ 0J + (a-l)$ , tJ $ ]^r 2 dr!. (22) 
Jov 4?r 

The first term is the rate of change of the external field 
energy ( |l8| ) inside V, resulting from the evolution of the 
tidal field. The second term is the rate of change of the 
interaction energy (pOl). The third and fourth terms to- 



gether, by comparison with Eq. (16), must be equal to 
dE /dt, the rate of change of the body's self energy. The 
fourth term is the contribution from the body's own field 
moving across the boundary dV. Therefore, the third 
term gives the change in the body's energy coming from 
the interaction with the tidal field; in other words, it is 
the rate of work done on the body by the tidal field. 
Furthermore, this term is independent of a and is, con- 
sequently, independent of how the Newtonian energy is 
localized, as claimed in Sec. I, Eq. (0). 



III. RELATIVISTIC ANALYSIS 

In this section, we will exhibit a relativistic version of 
the calculation in Sec. II, again showing that the rate 
of work done on the body by the tidal field is gauge 
invariant and that it has the same value in a general 
relativistic perturbative treatment as in the Newtonian 
one: — (l/2)£ijdTij/dt. The formalism to be used is 
the Landau-Lifshitz energy-momentum pseudotensor and 
multipole expansions as developed by Thorne and Hartle 
jjj and Zhang (SJ|(|, together with the slow- motion ap- 
proximation, so time derivatives are small compared to 
spatial gradients. 

We will work in the body's vacuum local asymptotic 
rest frame, which is defined as a region outside the body 
and far enough from it that its gravitational field can be 
regarded as a weak perturbation of flat spacetime, yet 
near enough that the tidal field of the external universe 
can be regarded as nearly homogeneous. This region is 
a spherical shell around the body; its inner boundary is 
near the body's surface but far enough away for gravity 
to be weak, and its outer boundary is at a distance where 
the external tidal field begins to depart from homogeneity 
(see Fig. 1). Somewhat more precisely, the local asymp- 
totic rest frame is the region throughout which r/C <C 1, 
r/lZ <C 1, and M/r <C 1, where r is the radial distance 
from the body, M is the mass of the body, and 1Z and C 
are the radius of curvature and the scale of homogeneity 
of the external gravitational field. If this region exists 
(as in the case of a black hole binary far from merger, 
for example) and the slow-motion limit applies, then the 
following analysis is valid. 

As in the Newtonian case, we will consider contribu- 
tions to the total energy inside a sphere V which encom- 
passes the body and whose boundary dV lies within the 
local asymptotic rest frame (see Fig. 1). 



A. deDonder Gauge 

We shall begin by computing the work done on the 
body by the external tidal field, using deDonder gauge 
(this section). Then in the next section, we shall show 
that the work is gauge invariant. 

DeDonder gauge is defined by the condition that 
1 ■- 0, where h a ^ is defined in terms of the metric 



density as follows: 
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(23) 



At linear order in the strength of gravity, h a P is the trace- 
reversed metric perturbation. According to Zhang ||, 
the components of h al3 in the body's local asymptotic 
rest frame are, at leading (linear) order in the strength of 
gravity and at leading non-zero order in our slow-motion 
expansion, 



h w = -4$ = 4— + 6 



(24) 



h^ = ~A 3 = 6^- + ^£ ab x a x b x> - ^r£ ja x a r 2 , (25) 



4 • 
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h ij = O 4-r 4 & =31 



(26) 



where the dots indicate time derivatives (i.e., = 
dlij/dt) and the symbol "&" means "plus terms of the 
form and magnitude." Note that the higher-order (£- 
order) multipoles have been omitted, since their con- 
tributions are smaller by ~ {r/C) 1 ^ 2 < r/C <C 1 and 
~ (M/r) l ~ 2 < M/r < 1 than the quadrupolar (£ = 2) 
terms that we have kept. The second time derivatives will 
also be neglected since they are unimportant in the slow- 
motion approximation; this effectively eliminates h lJ in 
this gauge. Also, note that the $ in Eq. ( p4] ) has the 
same form as the Newtonian gravitational potential of 
Sec. II, and the Aj of Eq. ( |25| ) is a gravitational vector 
potential, which does not appear in Newtonian theory. 

In general, it would be possible to have a term oc 
Tjk£jk/r in Eq. (pi|), as well as the cx M/r term. We have 
chosen to define the constant for the monopolar term to 
be 4M, thereby eliminating any term cx Xjk£jk/f] this is 
arbitrary but convenient, as will be seen shortly. 

The total mass-energy Aly inside the sphere V (to- 
tal stellar mass including the quadrupolar deformation 
energy and energy of interaction between the deforma- 
tion and tidal field) is defined by Thorne and Hartle (Eq. 
(2.2a) of Ref. 0) in terms of the Landau-Lifshitz super- 
potential as 



1 

16^ 
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where 



(27) 



(28) 
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cf. Eqs. (20.6), (20.3), and (20.20) of MTW {land Eqs. 
(100.14) and_(100.2) of Landau and Lifshitz |p) . Using 
Eqs. @, (§f), and (H)-(||) in Eq. (§?]) and carrying 
out the integration gives 



M v = M+O (sSr 1 & SSr 7 

+0 (±Sr 2 & ISr 2 & ISr 2 



„ / IT „ X J 
+0 & ^- 



(29) 



Hence, M.\> — M at leading order in our slow-motion 
approximation when we neglect the double and higher- 
order time derivatives. This is the reason the constant of 
the monopolar term in Eq. (|24|) was chosen to be AM. 

To calculate the rate of work done by the tidal field on 
the body when Sjk and Jjk change slowly, we consider 
the change of the mass-energy Mv, 



dMy 
dt 



(-g)t oj njr 2 dn 



(30) 



where t^ v is the Landau-Lifshitz pseudotensor. That this 
expression is indeed the time derivative of Eq. ( |27| ) is a 
result of Gauss's theorem (see discussion after Eqs. (2.3) 
ofRef. 0). 

In the body's local asymptotic rest frame, the Landau- 
Lifshitz energy-momentum pseudotensor (Eq. (20.22) in 
MTW H ) is, at the orders of accuracy we are considering, 

= , (31) 

= ^* A + 4^ ( A U ~ ■ (32) 



Here the deDonder gauge condition h g = nas been 
used to eliminate many terms from the general expression 
in MTW, but most of the simplification has come from 
keeping only terms of leading-order in the slow-motion 
approximation [zeroth and first time derivatives, respec- 
tively, in Eqs. ( pl| ) and (|32|)] . This restriction has given us 
only products of h fJlV a which will produce terms contain- 
ing the products M 2 , ME, Ml, IS, 11, SS for (-g)t 00 
and Ml, MS, IS, SI, IT, SS for (-g)t° 3 . This may 
be illustrated by expressing (— g)t 0a explicitly in terms of 
the quadrupole moments by substituting Eqs. ( |24| ) and 
(||) into Eqs. @ and © to get 



426 



Ml ab x a x b 
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^ab-^bc^ 



3 1 5 ^dfj^-QciX XXX 



(33) 



and 



/ \.fH 1 / -i n abScj% % 3* , ^.TabSbc'E % ^ 3 

(—g)t 3 = tzt- I -18 — h 24- 
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a^b^c q- .p. ^a^c^j 
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(34) 



Note that we have kept only the SI cross terms in the 
expression for {—g)t°^, as only they will contribute to our 
calculation of the interaction energy and work. 

We find the rate of change of mass-energy inside our 
sphere V by inserting Eq. (Q) into Eq. ([5(j|) and inte- 
grating. The result (only considering the cross terms) 



dM 



dt dt \ 10 13 13 



1„ d 

-f T. . 

2 13 dt %3 



(35) 



Since the interaction energy can depend only on the in- 
stantaneous deformation and tidal field, its rate of change 
must be a perfect differential, whereas the rate of work 
need not be. Also, if the tidal field changes but the body 
does not, there is no work done. From these two facts, we 
can conclude that the first term of Eq. (|3^) is the rate of 
change of the interaction energy between the tidal field 
and the body and that the second term represents the 
rate of work done by the external field on the body (the 
"tidal heating"). 

Notice that this value for the rate of work matches that 
discussed in Sec. I, Eq. (||), and found via the Newtonian 
analysis in Sec. II. Note the comparison with our New- 
tonian expressions (^2|) and (p0|). The first and fourth 
terms of Eq. ( p2| ) are missing here because we included 
in our (—g)t 0j only the (body field) x (external field) 
crossterms. If we had also included (external) x (external) 
and (body) x (body) terms, Eq. (H||) would have entailed 
expressions like the first and fourth terms of Eq. 
Note also from the interaction energy terms in Eqs. 
(pp|), and (35) that the Landau-Lifshitz way of localizing 
gravitational energy corresponds to the Newtonian choice 
a = — 3, a correspondence that has previously been de- 
rived by Chandrasekhar p5| . 
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B. General Gauge 

In Sec. III. A., we considered the special case of de- 
Donder gauge, which was particularly simple due to the 
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gauge condition and to the terms being effectively 
zero. Now we will examine a general gauge, which can 
be achieved by a gauge transformation of the form 



Sh/j,v <^i,i/+£i/,/j ~~ Vfiv£, ,a 



(36) 



where £ a is a function to be chosen shortly. Using 
h^v = h^v — ^TjfMuh^, where is the perturbation of 
the metric away from the Minkowski metric in the local 
asymptotic rest frame, this can also be expressed as 



Note that in the deDonder gauge, 



/ioo— 2- 



M 



Sij X X"^ , 



(37) 



(38) 



T T a 10 4 
S= - e£ 7^ - ^«6* a *V + ^£ 3a x a r 2 , (39) 



M l ke x k x e k i 

2 h3 = £ kt x x 



(40) 



Since we are interested only in gauge changes of the 
same order as we have been using so far (leading-order in 
the slow- motion approximation), we include only terms 
i n ( a that will produce Sh^ v of the same forms as Eqs. 
(^8|)-(^0|), but with different numerical coefficients. For 
example, consider Shoo — 2£o.o ^ M/r; that gives 
£o cx Mt/r, since M is a constant. This, in turn, im- 
plies Sh j — £o.j ^ Mtx 3 /r 3 , but this is not of the same 
form as the terms in Eq. (f39|); rather, it corresponds 
to a gauge that rapidly becomes ill-behaved as time 
passes. Similar arguments apply to Shoo oc Ij k x^ x k /r 5 
or Shoo 

OC £j X"^ X - their coefficients cannot be altered 
by a gauge change because such a change would alter the 
mathematical form of hoj and would make its magnitude 
unacceptably large in the slow-motion limit. As a result, 
we must set £o = 0. If we now consider Shoj = £j,o, pos- 
sible terms are of the form oc Ij a x a /r 3 , oc £ ab x a x b x J , or 
oc £ ja x a r 2 ; cf. Eq. (p9|). Each of these gives 5h Q o = 



and Shij of the same form as Eq. (40); hence, terms of 



these forms are allowed. Consequently, the most general 
gauge change that preserves the mathematical forms of 
h^ v but alters their numerical coefficients is 



r a 10 

h°i = (2 - a)X ja — + (— - 1 )£ ab x a x b x* 



(43) 



h jk = 2al jk — + SaS^lat—^- - 3al ja - 



-3al ka ^- + 2/3£ ]k r 2 - 2{f3 + 1 )5^ k £ ab x a x b 
+2((3 + 1 )£ ]a x a x k + 2(/3 + -f)£ ka x a xi . (44) 

In this general gauge, if we calculate the total mass- 
energy Mv inside the sphere V using Eqs. J2S|), (p3|), and 
@-@ in Eq. (§7]), we find 

/ 29 4 \ 

M v = M + f 2 7 2 + — /3 7 - - 7 + 2f3 2 - 2/3 J £ l3 £ l3 r 5 

1 7 4 23 \ 

-a + 2p + ~~/--ap-—aj)l ij £ ij . (45) 

The lijlij /r 5 term is zero. It is comforting to note that 
all the new terms vanish for a = (3 = 7 = 0, giving the 
deDonder result Mv = M . The ££ r 5 and T£ terms in 
Eq. ( |45| ) are large compared to the double time-derivative 
terms that formed the largest corrections to the mass- 
energy in deDonder gauge (|9|); however, they are still 
small compared to the mass M that appears in the ex- 
pansion (|^) °f h 00 . Also note that, near the body of 
interest, the ££r° term will be small compared to the I£ 
term, due to its radial dependence. So, once again, we 
have Mv ~ M as a first approximation, although it is 
necessary to keep the extra terms in Eq. (^) to main- 
tain the same level of accuracy as we have been using. 
Consequently, Mv is gauge-dependent to the order that 
interests us, and it has the "2jj£jj" ambiguity discussed 
by Thorne and Hartle [Q and mentioned in Sec. I. 

Keeping only the leading-order terms in the slow- 
motion approximation, as described in Sec. III. A., the 
Landau-Lifshitz pseudotensor in the new gauge is 



1 



7 



16tt \ 8 



.-h ij h lk ■ + -h Q0 h 3j 

o ,k >J 4 >* ,i 



k n ,k ' 



(46) 



6 = 0, 

& = a ^T~ + P£jax a r 2 + 1 £ ab x a x b x J , (41) 

where a, (3, and 7 are arbitrary constants. 

Using the trace-reversed gauge change (|6|) with Eq. 
([id]), the new h af3 become 



M 
r 



h m = 4— + (5 7 + 2/3- 2)£ ljX l x J 



+ (6 - 3a)l, 



v r 5 ' 



(42) 



+h°\Ji*\ l -h%h i \ j +h i \ k V\o 

*-ukk 100 1t00 Ifefc 



(47) 



7 



Note that the first term of Eq. (|46|) is the same as Eq. 
(|3l|), and the first three terms of Eq. ([47| ) are the same 
as Eq. (32). The additional terms all involve h? k , which 
was effectively zero in the dcDonder gauge because of our 
slow-motion assumption. 

Substituting Eqs. @ and ©-© into Eq. © and 
integrating gives the rate of change of mass-energy inside 
the sphere V as 



dt 



d 
~dl 



( 7aj 



97 



7a(3 



12/3 
~~5~ 



a 



1 

To 



1 dl.j 

2 ij ~dT 



(48) 



Again, we have kept only the (external field) x (body 
field) crossterms. As expected, this expression reduces 
to the deDonder result (jJTJ) when a = f3 = 7 = 0. Using 
the same argument for Eq. ( f48| ) as for Eq. ([35]), we can 
conclude that the first term of Eq. ([TJJ) is the rate of 
change of the interaction energy between the tidal field 
and the body and that the second term represents the 
rate of work done by the external field on the body (the 
"tidal heating"). 

Notice the gauge dependence (dependence on a, (3, 7) 
of the rate of change of interaction energy 



dE hl 
dt 



d 
~dt 



/7a-y 9 7 7a/3 12/3 a 1 
l~5 5~ + ~5 5~~ ~ 5 + 10 



(49) 



By contrast, the "tidal heating" work done on the body 
has the same, gauge-invariant value as in Newtonian the- 
ory 



dW 
~df 



1 „ dli 



— p 

- 2 1„ 



dt 



(50) 



IV. CONCLUSIONS 

In this paper, we have shown that the rate of work 
done by an external tidal field on a body is indepen- 
dent of how gravitational energy is localized in Newto- 
nian theory and that it is gauge invariant in general rela- 
tivity. Furthermore, this quantity — which we are calling 
the "tidal heating" — is given unambiguously by Eq. (||). 
That the tidal heating should be a well-defined and pre- 
cise quantity is reasonable, given that its physical effects 
have been observed in the form of volcanic activity of 
Jupiter's moon Io HJlJJlJ] . 

There remains one aspect of the uniqueness of the 
tidal heating that we have not explored: a conceivable 
(but highly unlikely) dependence of dWj dt on the choice 
of energy-momentum pseudotensor in general relativity. 
The arbitrariness of the pseudotensor is general relativ- 
ity's analog of Newtonian theory's arbitrariness of local- 
ization of gravitational energy. The fact that dW/dt is 



independent of the Newtonian energy localization sug- 
gests that it may also be independent of the general rel- 
ativistic pseudotensor. In addition, the clear physical 
nature of dW/dt gives further confidence that it must 
be independent of the pseudotensor. Nevertheless, it 
would be worthwhile to verify explicitly that dW/dt is 
pseudotensor- independent . 
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FIG. 1. An example of an isolated, slow-motion body: a 
star or black hole in a binary system with R/a <C 1, where R 
is the radius of the body and a is the separation of the body 
and companion. The dashed circles indicate the boundaries 
of the body's vacuum local asymptotic rest frame, the region 
in which M/r <C 1 and r/C <C 1. Here, r is the radial co- 
ordinate, M is the mass of the body, and £ is the scale of 
homogeneity of the gravitational field. The boundary of the 
sphere over which we integrate, denoted by dV, lies within 
this region. 
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